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INTRODUCTION 



. . is essentially a progress report on a theoretical 

investigation of the propagation of transient waves in a random 
medium. The emphasis in this study is on applications to 
sonic -ooom propagation, particularly as regards the effect of 
atmospheric turbulence on the sonic-boom waveform. The analysis 
is general, however, and is applicable to other types of waves 
besides sonic -boom waves . 

. . _ The Phenomenon of primary concern in this investigation is 
the f inestructure of the wave. The first figure shows what is 
meant by the f inestructure . 
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1 . PRELIMINARIES 


The uppermost sketch shows a typical transient waveform as 
it might look if observed under ideal conditions; i.e., in the 
absence of scattering, refraction, or any other mechanisms which 
might act to distort the waveform. This is what I call the 
incident wave. 

When a wave propagates in a real medium, such as the 
atmosphere, a random high-frequency structure, called the 
f inestructure, generally apears on the waveform, giving it the 
appearance shown in the bottom sketch. The f inestructure is 
defined formally, at least insofar as the present investigation 
is concerned, as the difference between the observed wave and its 
ensemble average, and generally has an appearance much like that 
shown in the middle sketch. 

The f inestructure is believed to arise as a result of 
scattering and redistribution of wave energy by the random 
inhomogeneities of the medium. The f inestructure is particularly 
important in connection with sonic-boom propagation, since the 
strong high-frequency component that often appears near the front 
of the N-wave may contribute appreciably to the perceived 
noisiness of the sonic boom. 

The f inestructure is most conveniently characterized by its 
statistical properties, such as the variance, the standard 
deviation, the correlation function, and the spectral density 
function. The spectral density function is the quantity that is 
most relevant to the question of annoyance, particularly as 
regards sonic-boom waves, since it contains information on the 
frequency content of the finestructure . The spectral density 
function is the quantity of primary interest in this study. 


< 


Incident Wave 





Observed Wave 



2 . FORMULATION 


The _ investigation is based on a relatively simple 
mathematical model. In this model the medium is assumed to be 
one-dimensional, with properties that vary only with the spatial 
coordinate; i.e., the medium is assumed to be independent of 
time. The medium is also assumed to be quiescent, which means 
that only thermal scattering is being considered. {Mechanical 
turbulence, which gives rise to inertial scattering, is ignored.) 
In addition, the medium is assumed to be non-dissipative . 
Non-linear effects are disregarded. 

The starting point of the analysis is the time -dependent 
scalar wave equation in one spatial dimension x. (Note that 
derivatives are denoted by superscripts rather than by 
subscripts. ) Here w is the wave function; f is the source term. 

The sound speed c*, which is assumed to be a random function 
of x, is written as shown, where c is a constant reference speed 
and fj, is a stationary random function of x having zero mean 
and unit variance. {The angle brackets denote an ensemble 
average.) _ The- paramter e is the standard deviation of the index 
of refraction of the medium, and is assumed to be small. 

It should be pointed out that, in the atmosphere, the 
standard deviation of the acoustic index of refraction is 
typically about one part in one thousand, so that the assumption 
that e is small, which is to say, the assumption that the medium 
is only weakly inhomogeneous, is appropriate for the study of 
sonic-boom propagation. 

„ , conditions reflect the assumption that the 

medium is at rest prior to the initial disturbance. 


c, 2 w tt -w xx = /; t > 0, -00 < X < + oo 

c 


c* (x) = 


1 — e/j,(x) 


, <fi>= 0 , < fjT >= 1 


w(x, 0) = w t (x,Q) = 0 
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3 . PERTURBATION METHOD 


The assumption of a weakly inhomogeneous medium allows the 
problem to be solved by a perturbation method. To begin, 
assume that the wave function w has an expansion (which may be 
asymptotic) in powers of c . By substituting this 
expansion back into the wave equation, equations for the 
coefficients wg , wj , etc. are derived in the usual way. The 
first two such equations, for wg and wi , are shown. It 
turns out to be sufficient to carry out this procedure only as 
far as the first -order term (the term) . The wave function w 
is then approximated by the first two terms of the expansion, as 
indicated on the last line. 

Instead of specifying the source term f, it turns out to be 
more convenient to specify the zeroth-order wave function; i.e., 
the Wg term, which corresponds to the incident wave. It is 
written as shown, where h, which is an arbitrary function, 
defines the waveform of the incident wave. 

Writing the incident wave in this form is equivalent 
to _ specifying a source term that is concentrated at the 
origin, having a time dependence determined by the function h. 

With Wg specified in terms of the function h, the right-hand 
side of the equation for wi is determined, provided that A 4 is 
regarded as a known function. That equation can be solved by the 
method of characteristics. That method is well known, and 
so will not be described here. Instead, let's 
look next at how an expression for the correlation 
function of the f inestructure is derived, once .an 
expression for wi has been obtained. 


w(x , t ; e) = wq(x , t) + ewi(x , t) 4- e 2 w^{x , t) + . . . 


2 n „tt „„xx r 

C W 0 - Wq = f 

—2 tt XX n —2 tt 

C Wi ~ Wf — 2 C {IWq 


wo(x,t) — h(t — \x\/c) 


w — it/g -f ew\ + 0(e 2 ) 
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4. CORRELATION FUNCTION 


_T° ^ et a general expression for the correlation function of 
the f inestructure , start with the expression for w and average 
it, noting that the average of w i is zero. (The average of wi is 
zero because it is linear in the random function fi , which has 
average zero. ) Note that the average field differs from the 
incident field by terms of order ^2 . 


Next, subtract to get an expression for w, the fluctuating 
rield, defined as the difference of w and its average. The 
quantity w corresponds to the f inestructure . An expression for 
the temporal correlation function of the f inestructure at any 
point x is then obtained by forming, at x, the product 
of w with itself at the two time values ti and t 2 • 

The averaging prcocedure used here is not the usual ensemble 
averaging • Instead, it's a travel— time-corrected averaging 
procedure that I call asynchronous ensemble averaging. The 
essential idea of the method is that, instead of measuring time 

r ® s P ect to some universal reference time, such as the time 
that the wave is emitted by the source, it is measured, for each 
wave in the ensemble, relative to the time that that wave arrives 
at the _ observer . The advantage of using a travel-time-corrected 
averaging procedure of this type is that it avoids certain 
spurious effects that arise as a consequence of averaging over an 
ensemble of waves that have become dispersed due to variations in 
travel time among the different members of the ensemble. 


A travel-time corrected averaging procedure similar to the 
one used here was used by Allan Pierce in one of his papers on 
sonic-boom propagation. 


£ should mention also that asynchronous averaging has the 
effect of renormalizing the perturbation series for the wave 
function, which results in the elimination of some secular terms 
from the expression for the first-order field. 


In order to apply asynchronous averaging in deriving an 
ex P^ ess ion for the correlation function of the f inestructure, one 
Pi k y t + r, where t is the travel time from the source 
(the origin)’ to the observation point x, and r is the time 
relative to the arrival time at which the observation is being 
made. The variable t 2 is replaced by t +■ r +■ s, where s denotes 
the separation time between the two observation times . The 
result is an expressicpn for the correlation function of the 
fines tructure , which is denoted by K. Note that terms of order 
have been dropped from the expression for K. 

The spectral density function is then obtained by taking the 
cosine transform of K with respect to s . 


( See next page . ) 
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CORRELATION FUNCTION 


w — wq + ewi 4- 0(e 2 ) 

< w > ~wq + 0(e 2 ) 
w— < w >= w = ewi 4- 0(e 2 ) 

< w(x,ti)w(x,t 2 ) > = e 2 < wi(x y ti)wi(x,t 2 ) > -K9(e 3 ) 

t\—r-\-r ; t 2 — r + r + s\ r > 0, s > 0 

K(x , r, 5 ) =< w(x, r 4- r)w(x , r 4- r 4- 5 ) > 

= e 2 < 4- r)wi(x,r 4- r* 4- 5 ) > 
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5. GENERAL RESULTS 


The result of the analysis sketched above is the general 
expression shown for the correlation function of the 
fines tructure . The function 7 that appears in this formula is a 
temporal correlation function that's related to the spatial 
correlation function of the medium by the change of variable 
indicated. The function h determines the waveform of the 
incident wave, as discussed previously; the primes denote 
derivatives. The lower limit of the second inner integral; i.e., 
f, is defined as shown. 

This formula is valid provided that the propagation path 
length x is much greater than the integral scale, or outer scale, 
L, of the medium. (This condition is satisfied in virtually all 
cases of practical interest.) Note also that this expression is 
valid only for r positive. When r is negative, i.e., before the 
wave has arrived, the correlation function is zero. 

forinula ' which expresses the correlation function of 
the fines tructure in terms of the incident waveform and the 
correlation functicpn of the medium (the last two quantities both 
being arbitrary), is the main result of this investigation. 

. The spectral density function of the f inestructure is 
obtained by taking the cosine transform of K with respect to s 
(the separation variable) . 

. One point about this expression that's worth noting is that 
it is independent of the propagation range x. The reason for 
^nis is that, according to the theory, by the time the incident 
wave has propagated a distance of the order of an outer scale 
length into the medium the f inestructure has become, in a 
statistical sense, fully developed. Beyond this range, i.e., for 
x >> L, there is, for all practical purposes, no further 
evolution of the statistical properties of the f inestructure . 

The correlation function K is, however, strongly dependent 
on the magnitude of the randomness of the medium, as is shown by 
the presence of the term e 2 on the right-hand side of this 
expression. 

Note also _ that the correlation function K depends on r (the 
elapsed time since onset) , as well as on s, the separation, time. 
The spectral density function of the f inestructure will therefore 
a ^ so on r * The correlation and spectral density functions 

or the f inestructure are thus time-dependent functions . This is 
a consequence of the fact that the finestructure of a transient 
signal is generally a non-stationary random process . 

The formula shown expresses the function K as the sum of two 
iterated integrals . Of these, the two inner integrals involve 
only the function h. If that function has a sufficiently simple 
form, then the inner integrals can be evaluated analytically, 
after which only a pair of single integrals is left to deal with. 

One case in which this type of simplification is feasible is 
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when the incident wave has the form of a simple ramp function, 
since then the function h is piecewise constant. The relevant 
calculations have been carried out, and the results are shown on 
the next figure. (This type of simplification is also feasible 
when the incident wave has the form of an N-wave, but the 
calculations for that case have not been done.) 


GENERAL RESULTS 


K(r y s) 



h* (77 - f)/i' (77 -4- s)d rjd^ 
h' (77 - g)h' (r))dr)dt\ 


7(0 = < H S )H S + 0 >1 <P(s) = fi Q 
£ = max( 0 , f), x » L, r > 0 , s > 0 
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6 . RAMP FUNCTION 


The sketch shows the function h(t) for the case of a 
ramp- function incident wave. The amplitude P and the rise time 6 
are shown on the sketch, as are the two time values r and r+s at 
which the fmestructure is evaluated in order to form the 
two -point correlation function. The results described here are 
for the case in which the points r and r+s are both behind the 
rise phase of the incident wave, as indicated. This is the case 
of most practical importance. 


As has been mentioned, the two inner integrals in the 
expression for K shown on the previous vu-graph can, in this 
case, be evaluated exactly. The result of that calculation is 
the first equation shown. Note that, although K is written 
generally as a function of both r and s, the right-hand side of 
this expression is in fact independent of r, showing that the 
fmestructure is a stationary stochastic process in this case. 

An expression for the spectral density function of the 
fmestructure, denoted by D, is obtained by taking the cosine 
transform of K with respect to s. The result is given by the 
second equation. The function is the transform of 7 . 

Of these two quantities, the spectral density function ,- 
i.e., D, t is, as has already been emphasized, the quantity of 
primary interest in this study. The expression for this quantity 
will be examined more closely in a moment. First, however, 
let s ( take a quick look at what happens to these two formulas in 

^i e as S goes to zero, i.e., as the ramp function becomes a 

step function. 


The limiting forms of K and D, denoted by K 0 (s) and D 0 (a/) , 
respectively, are easy to calculate, and are given by the last 
two equations. The function Kq(s) is then the correlation 
function, while Dq {u) is the spectral density function, of the 
fmestructure for the case in which the incident wave is a step 
function. 

It should be pointed out that no analytical problems arise 
in calculating _ these limits. There is thus no problem with 
singular behavior of the f inestructure in the limit as the ramp- 
tunction incident wave becomes a step function. 

These last two formulas show that, when the incident wave is 
a step function, the correlation function of the f inestructure is 
proportional to the correlation function of the medium, and, 
similarly , the spectral density function of the f inestructure is 
proportional to the spectral density function of the medium. 

When the incident wave has the form of a ramp function, 
however, with a non-zero rise time, these simple relations no 
longer obtain. This situation, particularly as regards the 
spectral density function, will be examined in more detail in 
the next section. 
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7. RAMP FUNCTION (cont'd) 


t In order to examine some of the implications of the result 
(Detained above for the ramp-function f inestructure spectrum D, 
it's convenient to write the expression for that quantity in the 
form shown on the first line. In this formula the function Dn is 
the f inestructure spectrum for the case of zero rise time, i.e., 
for the step-function incident wave, as described above. The 
function A determines how the f inestructure spectrum differs from 
Dq when the rise time 6 is non-zero . 

The function A is sketched in the figure. Note that it is 
effectively zero when its argument is greater than about 27T. What 
this means is that (referring to the formula for D) the 
fines trmeture spectrum is effectively zero at all frequencies for 
which glu is greater than about 27 r, or, what is equivalent, when 
f, the frequency in Hz, is greater than about 1/6 . 

. This result shows that, the greater the rise time of. the 
incident wave, the less high-frequency energy there is in the 
f inestructure . A couple of examples of this effect, involving 
parameter values typical of sonic-boom propagation, are shown: 

When the rise time is one millisecond, there is effectively no 
energy in the f inestructure spectrum at frequencies above about 1 
kHz. If the rise time is increased to two milliseconds, then the 
lower limit of the frequency range for which the f inestructure 
spectrum is devoid of energy drops to 500 Hz. 

Note, however, that, since the function A is equal to unity 
w i3 en 4 r : ^ s argument is equal to zero, the low-frequency portion of 
the f inestructure spectrum is relatively unaffected by changes in 
the rise time. 

Since it's generally the high-frequency portion of the 
acoustic spectrum that's most annoying to the human ear, what 
these > results imply is that, at least for a ramp-function wave in 
one dimension, increasing the rise time of the incident wave will 
tend to reduce the annoyance associated with the f inestructure . 
Whether this result holds in three dimensions as well remains to 
be seen, but if it does it has obvious implications as regards 
the idea of shaping the sonic boom in order to reduce its 
annoyance . 

Th® results just described could, of course, be expressed 
just as well in terms of time scales, rather than in terms of 
frequency. _ Expressed in those terms, the results imply that a 
ramp- function wave propagating in a random medium will contain 
very little structure having time scales smaller than the rise 
time. This prediction agrees with observations of sonic-boom and 
other types of transient waves, which rarely show any appreciable 
structure having time scales smaller than the rise time. 
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RAMP FUNCTION (cont'd) 


D{r,w) = A(Soj) D a {u) 

The function A(8u) r and therefore the f inestructure spectrum, 
effectively zero for 8u > 2tt , or f>l/6. 

Sonic boom: 8 = .001 sec; /> 1000 Hz 

8 = .002 sec; /> 500 Hz 
D(r,uj) ~ 0 


8 . SUMMARY 


fnlln The ^ain results of this investigation can be summarized as 
f lrst ' ^he approach described here has been found to 
G ? nC * Se ' general expressions for the correlation and 
spectral density functions of the f inestructure . Second, the 
results for the case of a ramp- function incident wave agree 

S ? bservatlon that the important time scales associated 
with the f inestructure of sonic-boom and other types of transient 

Finall5 r %h enera1 }^ GG ™Parable to, or greater than, the rise time. 

^t^ reSUltS (again fo 5 the ramp-function incident wave) 

• t !3 a 5 an ?- nGrease 111 the rise time of the incident wave 

the f inestructure a reductlon in the high-frequency content of 

T his last result, assuming that it holds in the 
three-dimensional case as well, has implications as regards the 
iaea ot shaping the sonic boom in order to reduce its annoyance 

? i? eSta 5 hat incr easing the rise time of the sonic boom will 
maxe it quieter, even m the presence of turbulence. 
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